Association of genetic vatiation within UBL5 with phenotypes of metabolic syndrome

The BEACON gene was initially identified using the differential display polymerase chain reaction on hypothalamic mRNA samples collected from lean and obese Psammomys obesus, a polygenic animal model of obesity. Hypothalamic BEACON gene expression was positively correlated with percentage of body fat, and intracerebroventricular infusion of the Beacon protein resulted in a dose-dependent increase in food intake and body weight. The human homolog of BEACON, UBL5, is located on chromosome 19p in a region previously linked to quantitative traits related to obesity. Our previous studies showed a statistically significant association between UBL5 sequence variation and several obesity- and diabetes-related quantitative physiological measures in Asian Indian and Micronesian cohorts. Here we undertake a replication study in a Mexican American cohort where the original linkage signal was first detected. We exhaustively resequenced the complete gene plus the putative promoter region for genetic variation in 55 individuals and identified five single nucleotide polymorphisms (SNPs), one of which was novel. These SNPs were genotyped in a Mexican American cohort of 900 individuals from 40 families. Using a quantitative trait linkage disequilibrium test, we found significant associations between UBL5 genetic variants and waist-to-hip ratio (p = 0.027), and the circulating concentrations of insulin (p = 0.018) and total cholesterol (p = 0.023) in fasted individuals. These data are consistent with our earlier published studies and further support a functional role for the UBL5 gene in influencing physiological traits that underpin the development of metabolic syndrome.

Development of twenty-four novel microsatellite markers for the freshwater crayfish, Geocharax gracilis, using next generation sequencing

The crayfish Geocharax gracilis is an important inhabitant of natural and agricultural drainage systems in south-eastern Australia. To investigate population structure, genetic diversity and patterns of connectivity in natural and human-altered ecosystems, we isolated and characterised 24 microsatellite loci using next generation sequencing. Loci were initially tested for levels of variation based on 12 individuals from across the species’ geographical range. A further 33 individuals from a single population were used to test for departures from Hardy–Weinberg equilibrium and linkage disequilibrium. We detected high to moderate levels of genetic variation across most loci with a mean allelic richness of 8.42 and observed heterozygosity of 0.629 (all samples combined). We found no evidence for linkage disequilibrium between any loci and only three loci (Geo01, Geo24 and Geo47) showed significant deviations from Hardy–Weinberg expectations. These same three loci, plus two additional loci (Geo06 and Geo28), also showed the presence of null alleles. These 24 variable markers will provide an important tool for future population genetic assessments in natural and human altered environments.

Development of twenty-three novel microsatellite markers for the seagrass, Zostera muelleri from Australia

Seagrasses are one of the most productive and economically important habitats in the coastal zone, but they are disappearing at an alarming rate, with more than half the world’s seagrass area lost since the 1990s. They now face serious threat from climate change, and there is much current speculation over whether they will survive the coming decades. The future of seagrasses depends on their ability to recover and adapt to environmental change—i.e. their ‘resilience’. Key to this, is understanding the role that genetic diversity plays in the resilience of this highly clonal group of species. To investigate population structure, genetic diversity, mating system (sexual versus asexual reproduction) and patterns of connectivity, we isolated and characterised 23 microsatellite loci using next generation sequencing for the Australian seagrass species, Zostera muelleri (syn. Z. capricorni), which is regarded as a globally significant congeneric species. Loci were tested for levels of variation based on eight individuals sampled from Lake Macquarie, New South Wales, Australia. We detected high to moderate levels of genetic variation across most loci with a mean allelic richness of 3.64 and unbiased expected hetrozygosity of 0.562. We found no evidence for linkage disequilibrium between any loci and only three loci (ZosNSW25, ZosNSW2, and ZosNSW47) showed significant deviations from Hardy–Weinberg expectations. All individuals displayed a unique multi-locus genotype and the combined probability of identity across all loci was low (P _ID = 1.87 × 10⁻¹²) indicating a high level of power in detecting unique genotypes. These 23 markers will provide an important tool for future population genetic assessments in this important keystone species.

The development of 10 novel polymorphic microsatellite markers through next generation sequencing and a preliminary population genetic analysis for the endangered Glenelg spiny crayfish, Euastacus bispinosus

The Glenelg spiny crayfish, Euastacus bispinosus, is an iconic freshwater invertebrate of south eastern Australia and listed as 'endangered' under the Environment Protection and Biodiversity Conservation Act 1999, and 'vulnerable' under the International Union for Conservation of Nature's Red List. The species has suffered major population declines as a result of over-fishing, low environmental flows, the introduction of invasive fish species and habitat degradation. In order to develop an effective conservation strategy, patterns of gene flow, genetic structure and genetic diversity across the species distribution need to be clearly understood. In this study we develop a suite of polymorphic microsatellite markers by next generation sequencing. A total of 15 polymorphic loci were identified and 10 characterized using 22 individuals from the lower Glenelg River. We observed low to moderate genetic variation across most loci (mean number of alleles per locus = 2.80; mean expected heterozygosity = 0.36) with no evidence of individual loci deviating significantly from Hardy-Weinberg equilibrium. Marker independence was confirmed with tests for linkage disequilibrium, and analyses indicated no evidence of null alleles across loci. Individuals from two additional sites (Crawford River, Victoria; Ewens Ponds Conservation Park, South Australia) were genotyped at all 10 loci and a preliminary investigation of genetic diversity and population structure was undertaken. Analyses indicate high levels of genetic differentiation among sample locations (F ST = 0.49), while the Ewens Ponds population is genetically homogeneous, indicating a likely small founder group and ongoing inbreeding. Management actions will be needed to restore genetic diversity in this and possibly other at risk populations. These markers will provide a valuable resource for future population genetic assessments so that an effective framework can be developed for implementing conservation strategies for E. bispinosus.

Microsatellite loci and the complete mitochondrial DNA sequence characterized through next generation sequencing and de novo genome assembly for the critically endangered orange-bellied parrot, Neophema chrysogaster

A suite of polymorphic microsatellite markers and the complete mitochondrial genome sequence was developed by next generation sequencing (NGS) for the critically endangered orange-bellied parrot, Neophema chrysogaster. A total of 14 polymorphic loci were identified and characterized using DNA extractions representing 40 individuals from Melaleuca, Tasmania, sampled in 2002. We observed moderate genetic variation across most loci (mean number of alleles per locus = 2.79; mean expected heterozygosity = 0.53) with no evidence of individual loci deviating significantly from Hardy-Weinberg equilibrium. Marker independence was confirmed with tests for linkage disequilibrium, and analyses indicated no evidence of null alleles across loci. De novo and reference-based genome assemblies performed using MIRA were used to assemble the N. chrysogaster mitochondrial genome sequence with mean coverage of 116-fold (range 89 to 142-fold). The mitochondrial genome consists of 18,034 base pairs, and a typical metazoan mitochondrial gene content consisting of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a single large non-coding region (control region). The arrangement of mitochondrial genes is also typical of Avian taxa. The annotation of the mitochondrial genome and the characterization of 14 microsatellite markers provide a valuable resource for future genetic monitoring of wild and captive N. chrysogaster populations. As found previously, NGS provides a rapid, low cost and reliable method for polymorphic nuclear genetic marker development and determining complete mitochondrial genome sequences when only a fraction of a genome is sequenced.

Characterization of microsatellite DNA markers for a mahseer species, Tor tambroides (Cyprinidae) and cross-amplification in four congeners

This study reports the isolation and characterization of microsatellite DNA markers in a mahseer species, Tor tambroides (Pisces, Cyprinidae). Of a total of 14 loci evaluated, 10 were polymorphic in T. tambroides samples, with an average of 2.86 alleles per locus. Deviations from Hardy–Weinberg equilibrium were observed at one locus and there was no indication of linkage disequilibrium among loci. A high level of cross-amplification among four congeners was achieved, with 12 loci successfully amplifying and 11 loci showing polymorphism in at least one other species. These markers will be a useful resource for population genetic studies and broodstock management of closely related mahseer species.

Antihypertensive treatments obscure familial contributions to blood pressure variation

The linkage and association between inherent blood pressure and underlying genotype is potentially confounded by antihypertensive treatment. We estimated blood pressure variance components (genetic, shared environmental, individual-specific) in 767 adult volunteer families by using a variety of approaches to adjusting blood pressure of the 244 subjects (8.2%) receiving antihypertensive medications. The additive genetic component of variance for systolic pressure was 73.9 mm Hg(2) (SE, 8.8) when measured pressures (adjusted for age by gender within each generation) were used but fell to 61.4 mm Hg(2) (SE, 8.0) when treated subjects were excluded. When the relevant 95th percentile values were substituted for treated systolic pressures, the additive genetic component was 81.9 mm Hg(2) (SE, 9.5), but individual adjustments in systolic pressure ranged from -53.5 mm Hg to +64.5 mm Hg (mean, +17.2 mm Hg). Instead, when 10 mm Hg was added to treated systolic pressure, the additive genetic component rose to 86.6 mm Hg(2) (SE, 10.1). Similar changes were seen in the shared environment component of variance for systolic pressure and for the combined genetic and shared environmental (ie, familial) components of diastolic pressure. There was little change in the individual-specific variance component across any of the methods. Therefore, treated subjects contribute important information to the familial components of blood pressure variance. This information is lost if treated subjects are excluded and obscured by treatment effects if unadjusted measured pressures are used. Adding back an appropriate increment of pressure restores familial components, more closely reflects the pretreatment values, and should increase the power of genomic linkage and linkage disequilibrium analyses.

The development of 20 microsatellite loci for the Australian marine mollusk, Donax deltoides, through next generation DNA sequencing

 Next Generation DNA sequencing was used to develop a suite of microsatellite markers for the marine mollusk, Donax deltoides. A total of 20 polymorphic loci were identified and 12 characterized using 30 individuals from a single population (Venus Bay) in south eastern Australia. We observed moderate to high genetic variation across most loci (mean number of alleles per locus = 7.3; mean heterozygosity = 0.633) with only a single locus (Ddel32) displaying significant deviation from Hardy–Weinberg equilibrium. Marker independence was confirmed with tests for linkage disequilibrium, however two loci were found to be influenced by null alleles. The 10 viable markers characterized in the present study provide a valuable resource for future population genetic assessments and fisheries management of D. deltoides in Australia.

ABCB1 polymorphism predicts escitalopram dose needed for remission in major depression

The ATP-binding cassette family of transporter proteins, subfamily B (MDR/TAP), member 1 (ABCB1) (P-glycoprotein) transporter is a key component of the blood–brain barrier. Many antidepressants are subject to ABCB1 efflux. Functional polymorphisms of ABCB1 may influence central nervous system bioavailability of antidepressants subject to efflux. Single-nucleotide polymorphisms (SNPs) at rs1045642 (C3435T) of ABCB1 have been associated with efflux pump efficiency. This may explain part of the interindividual variation in antidepressant dose needed to remit. Individuals (N=113) with DSM-IV (Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition) major depressive disorder (MDD) were treated with escitalopram (ESC) or venlafaxine (VEN) over 8 weeks. The17-item Hamilton Depression Rating Scale was assessed serially, blind to genotype. SNP rs1045642 of ABCB1 along with two SNPs previously reported to be in linkage disequilibrium with it (rs2032582 and rs1128503) were genotyped. Demographic features, clinical features, P450 metabolizer status and 5-HTTLPR (serotonin-transporter-linked promoter region) genotype were controlled for. Carriers of rs1045642 TT needed on average 11 mg of ESC to remit, whereas TC and CC carriers required 24 and 19 mg, respectively (P=0.0001). This equates to a 2.0- (95% confidence interval=1.5–3.4; P<0.001) fold greater ESC dose needed to remit for C carriers compared with TT carriers at rs1045642. Of VEN-treated subjects carrying TT genotype at rs1045642, 73.3% remitted compared with 12.5% for CC genotype (odds ratio=6.69; 95% confidence interval=1.72–25.9, P=0.006). These data suggest that antidepressant dose needed to remit can be predicted by an ABCB1 SNP. This has the potential clinical translation implications for dose selection and remission from MDD.

Characterization of nine polymorphic microsatellite loci in the dyeing poison frog Dendrobates tinctorius (Dendrobatidae), and their cross-species utility in two other dendrobatoid species

While field and laboratory based studies have provided significant insights into the parental care and courtship behaviour of dendrobatoid frogs, a comprehensive assessment of their genetic mating systems and population genetic parameters has been precluded because ofthe lack of highly variable DNA markers. Here we document the development of nine novel polymorphic microsatellite markers for the dyeing poison frog Dendrobates tinct or ius (Dendrobatidae ). We found between three and 16 alleles per locus in 60 individuals (30 males, 30 females) from the field site Saut Parare, French Guiana, with an average observed heterozygosity of 0. 75. None of the loci deviated significantly from Hardy-Weinberg equilibrium or showed linkage disequilibrium. We also report successful cross-species amplification of the nine markers in two other dendrobatoid species (Allobates femora/is and Oophaga pumilio). These markers have the potential to aid in determining the genetic structure of local populations, identifying small-scale phylogenies such as parent-offspring relationships and will allow for cross-species comparisons within dendrobatoid species. Therefore, these markers can be applied to a wide range of scientific fields, such as conservation, behavioural ecology and evolutionary biology.

Characterisation of the complete mitochondrial genome and 13 microsatellite loci through next-generation sequencing for the New Caledonian spider-ant Leptomyrmex pallens

The complete mitochondrial genome and a set of polymorphic microsatellite markers were identified by 454 pyrosequencing (1/16th of a plate) for the New Caledonian rainforest spider-ant Leptomyrmex pallens. De novo genome assembly recovered the entire mitochondrial genome with mean coverage of 8.9-fold (range 1-27). The mitogenome consists of 15,591 base pairs including 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a non-coding AT-rich region. The genome arrangement is typical of insect taxa and very similar to the only other published ant mitogenome from the Solenopsis genus, with the main differences consisting of translocations and inversions of tRNAs. A total of 13 polymorphic loci were also characterized using 41 individuals from a single population in the Aoupinié region, corresponding to workers from 21 nests and 16 foraging workers. We observed moderate genetic variation across most loci (mean number of alleles per locus = 4.50; mean expected heterozygosity = 0.53) with evidence of only two loci deviating significantly from Hardy-Weinberg equilibrium due to null alleles. Marker independence was confirmed with tests for linkage disequilibrium. Most loci cross amplified for three additional Leptomyrmex species. The annotation of the mitogenome and characterization of microsatellite markers will provide useful tools for assessing the colony structure, population genetic patterns, and dispersal strategy of L. pallens in the context of rainforest fragmentation in New Caledonia. Furthermore, this paper confirms a recent line of evidence that comprehensive mitochondrial data can be obtained relatively easily from small next-generation sequencing analyses. Greater synthesis of next-generation sequencing data will play a significant role in expanding the taxonomic representation of mitochondrial genome sequences.

Obesity and diabetes gene discovery approaches

New treatments are currently required for the common metabolic diseases obesity and type 2 diabetes. The identification of physiological and  biochemical factors that underlie the metabolic disturbances observed in obesity and type 2 diabetes is a key step in developing better therapeutic outcomes. The discovery of new genes and pathways involved in the  pathogenesis of these diseases is critical to this process, however  identification of genes that contribute to the risk of developing these diseases represents a significant challenge as obesity and type 2 diabetes are complex diseases with many genetic and environmental causes. A number of diverse approaches have been used to discover and validate potential new targets for obesity and diabetes. To date, DNA-based approaches using candidate gene and genome-wide linkage analysis have had limited success in identifying genomic regions or genes involved in the development of these diseases. Recent advances in the ability to evaluate linkage analysis data from large family pedigrees using variance components based linkage analysis show great promise in robustly identifying genomic regions associated with the development of obesity and diabetes. RNA-based technologies such as cDNA microarrays have identified many genes differentially expressed in tissues of healthy and diseased subjects. Using a combined approach, we are endeavouring to focus attention on differentially expressed genes located in chromosomal regions previously linked with obesity and / or diabetes. Using this strategy, we have identified Beacon as a potential new target for obesity and diabetes.

VNTR polymorphism of the insulin gene and childhood overweight in a general population

Objective: The VNTR polymorphism 5' of the insulin gene has been related to obesity in a previous study on children with early onset of severe obesity. Our purpose was to analyze the association between this polymorphism and adiposity variability in an unselected population of children and adolescents in northern France.

Research Methods and Procedures: In 293 nuclear families from the Fleurbaix Laventie Ville Santé study, we genotyped the INS VNTR polymorphism in 431 children and adolescents (8 to 18 years of age) and their parents. Overweight was defined according to the international definition in both children and adults. A transmission disequilibrium test in families with an overweight offspring was performed. The prevalence of overweight was compared according to genotype. The effect of the genotype on BMI and waist circumference was tested with a linear regression model, adjusting for age, gender, and Tanner stage.

Results: There was an undertransmission of class III alleles from heterozygous parents to their overweight offspring (p < 0.002). Overweight was associated with class I alleles in children and adolescents (12% I/I, I/III vs. 3% III/III; p < 0.08). Those with a class III/III genotype had a 1 kg/m2 lower mean BMI (p = 0.04) and 3 cm lower waist circumference (p = 0.02) than those bearing one or two class I alleles. No association of adiposity or obesity with class I alleles was found in parents.

Discussion: INS VNTR polymorphism seems to contribute to differences in adiposity level in the general population of children and adolescents.

Analysis of protein sequence and interaction data for candidate disease gene prediction

Linkage analysis is a successful procedure to associate diseases with specific genomic regions. These regions are often large, containing hundreds of genes, which make experimental methods employed to identify the disease gene arduous and expensive. We present two methods to prioritize candidates for further experimental study: Common Pathway Scanning (CPS) and Common Module Profiling (CMP). CPS is based on the assumption that common phenotypes are associated with dysfunction in proteins that participate in the same complex or pathway. CPS applies network data derived from protein–protein interaction (PPI) and pathway databases to identify relationships between genes. CMP identifies likely candidates using a domain-dependent sequence similarity approach, based on the hypothesis that disruption of genes of similar function will lead to the same phenotype. Both algorithms use two forms of input data: known disease genes or multiple disease loci. When using known disease genes as input, our combined methods have a sensitivity of 0.52 and a specificity of 0.97 and reduce the candidate list by 13-fold. Using multiple loci, our methods successfully identify disease genes for all benchmark diseases with a sensitivity of 0.84 and a specificity of 0.63. Our combined approach prioritizes good candidates and will accelerate the disease gene discovery process.

Conservation of the ST6Gal I gene and its expression in the mammary gland

Milk sialoglycoconjugates can protect the gastrointestinal tract of the suckling neonate by competitively binding to invading pathogens and promoting growth of beneficial flora, and their potential role in postnatal brain development is of particular interest in human infant nutrition. Although the concentration and the distribution of sialoglycoconjugates have been extensively studied in the milk of various species, the investigation of sialyltransferase gene expression in the mammary gland, in the context of lactation, has been limited. The sialyltransferase enzyme ST6Gal I transfers sialic acid from CMP-sialic acid to type 2 (Galβ1,4GlcNAc) free disaccharides or the termini of N- or O-linked oligosaccharides using an α2,6-linkage. Expression of the ST6Gal I gene is primarily regulated at the level of transcription through the use of several cell and development- specific promoters, producing transcripts with divergent 5′ untranslated regions (UTR). In the mouse mammary gland, the novel 5′UTR exon (L) appears to be associated with a drastic increase in ST6Gal I gene expression during lactation. We find that rats also possess an exon (L), suggesting conservation of this regulatory mechanism in rodents. In contrast, an exon (L)-containing transcript was not detected in the lactating bovine or human mammary gland. We also observed a trend of increasing ST6Gal I gene expression in the bovine mammary gland, culminating in involution. This is in contrast to species such as mice where the greatest change in ST6Gal I gene expression occurs between pregnancy and lactation, suggesting different roles in rodents vs. other mammals for α2,6-sialylated oligosaccharides present in milk.